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value times an exponential tail. 
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1 Introduction 



Random matrix theory has important apphcations quite generally in the sta- 

tistical description of complex systems as e.g. for complex nuclei, for which it has 
been originally developed, or for chaotic systems with just a few degrees of freedom 
as treated in quantum chaos. Usually one restricts oneself to Gaussian ensembles of 
random matrices, meaning that the matrix elements have a Gaussian distribution 
where the diagonal matrix elements have all the same dispersion, the off-diagonal 
elements also have the same dispersion, but the former is by a factor 2 larger than 
the latter one. See the remarks in section 2. This ensemble is the only one which 
is invariant under symmetry transformations of the underlying matrix A. If A is 
real and symmetric, the group of relevant transformations consists of the orthog- 
onal transformations and we speak of the Gaussian Orthogonal Ensemble (GOE), 
while for Hermitian A the relevant group of transformations is that of the unitary 
transformations and we speak of the Gaussian Unitary Ensemble (GUE) of random 
matrices. The property of the matrix element distribution to be Gaussian is a di- 
rect consequence of just two assumptions, namely statistical independence of the 
distributions of the matrix elements and invariance against the group of appropriate 
transformations. We usually have in mind infinite dimensional matrices, although 
for many purposes finite dimensionality is useful and sufficient. 

Several generalizations are possible. One is the generalization towards general non- 
normal matrices, either fully complex [6] (see also [1]) but still Gaussian (invariant), 
or real positive but no longer invariant against the above group of transformations 
while still Gaussian (with different variances for different matrix elements) [7]. In 
the latter paper [7] we have also treated 2D real symmetric matrices whose matrix 
elements are Gaussian distributed but with different variances of the diagonal and 
the non-diagonal matrix elements. Thus they no longer enjoy GOE invariance, but 
still have Gaussian distributed matrix elements. 

In the present paper we study the level spacing distribution P{S) of 2D real random 
matrices by considering general distributions of the matrix elements, going beyond 
Gaussian. Therefore these ensembles of random matrices are no longer invariant 
under the mentioned transformation groups. The statistics P{S) of the level spacings 
S changes under transformations with the elements of those groups, so it depends 
on the basis chosen for their representation. Nevertheless, they are important in 
specific physical situations and also involve interesting mathematics. 

We shall first treat rigorously the case of general 2D non-normal matrice^ with 
Gaussian distributed matrix elements even admitting different variances of the vari- 
ous matrix elements. The level repulsion exponent p in this case will turn out to be 
still p = 1, unless there are constraints on the matrix elements like e. g. having only 

^ A matrix A is non-normal [T5] if it does not commute with its adjoint, i. e., [A, A+] ^ 0. 
Non- normal matrices have important physical applications, especially in dissipative systems [7], 

m, m, m, m- 
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positive or only negative nondiagonal elements etc. More general matrix element 
distributions in this case are left for future research, as they introduce additional 
serious mathematical problems, cf. section [21 

Second we consider symmetric matrices with other, general distributions of the real 
matrix elements. For such symmetric, real matrices we shall show that the level 
repulsion exponent is always p = 1, provided the distributions of all matrix elements 
are regular at zero value. Next we study in detail and without approximations the 
following specific cases of such regular distributions, namely: the uniform (box) 
distribution, the Cauchy-Lorentz distribution, the exponential distribution, and also 
some addenda to the Gaussian case, already dealt with in [7j. 

If, in contrast, the distribution of matrix elements is singular at zero value, P{S) 
shows different behaviour at S* = 0. For example, if the singularity of the statistics 
of the matrix elements is an integrable power law at zero value of S, the level 
spacing statistics P{S) exhibits a fractional exponent pawer law level repulsion as 
was discovered and treated in [H],[S], [ID], [II]- This probably is characteristic for 
sparsed matrices, which in turn are important random matrix models for nearly 
integrable systems of the KAM type. 



2 General real 2D matrices 



2.1 Level spacing statistics 



Consider 2x2 matrices A = [Aij), where i,j = l or 2. This matrix has two diagonal 
elements, which can always be chosen as a and —a. For general Aij introduce 
As = ^ {All + A22) and subtract the diagonal matrix Agl. Then An — Ag = a = 
— {A22 — As), i. e., one obtains the formula ([T]) without loss of generality. Quite 

generally, for a matrix ^ ^ ^ ^ the level spacing S = |Ai — A2I = |-y/(a — dy + 46c | 

only depends on the difference a — d, so that we can arbitrarily shift a and d by a 
constant, in particular by A^. H 

Let a as well as the nondiagonal elements bi and 62 all be real and write 



A ^ (A,) - ( :^ \ ) . (1) 

If 61 = 62, the matrix A is symmetric. 



^Let us make clear that in the general symmetric GOE matrix ^ ^ J the variances of the 

diagonal elements a and d are equal, but by a factor 2 larger than the variance of the offdiagonal 
element b. However, setting d = —a implies, that the GOE case occurs when the variance of a is 
equal to the variance of h. See also subsection 4.1 
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The eigenvalues of A follow from 
|A- All = 



a — A 61 
62 —a — A 



A^ - - 6162 = 



(2) 



i.e., Ai,2 = ±\/a^ + &1&2 • (3) 

The eigenvalues are real for arbitrary symmetric (61 = 62) matrices. In the more 
general case of hi 7^ 62 the eigenvalues are still real, if the product 6162 is larger than 
—a? , otherwise they are purely imaginary, but never general complex. 

If the matrix A is not symmetric, it is no longer normal. Namely, in the general 
case one finds for the commutator 



[A,A^ 



ft2 - /)2 



20(62 — bi 



2a{b2 - bi) bl - bl 



(4) 



Apparently the commutator [A, A+] = is zero iff 62 = 61, i- e., in the symmetric 
case. In general, 2x2 matrices A are non-normal, [A, A+] 7^ 0. 

The distribution P{S) of the level spacings 5* is given by 



/+00 c+00 r+oo / 
/ / dadbi db2 5[S -2 
-oo J —00 v/— 00 \ 



+ 6162 



9a{a)gbAbi)9b2{b2)- (5) 



Here 5{.) is the Dirac delta function and ga{o), Qhiibi), gh2{b2) are the normalized 
probability density functions for the matrix elements 0,61,62, respectively. P{S) is 
the central object of our study in this paper. We are going to study the dependence 
of the main properties of P{S), especially the small-S" behaviour (the level repulsion) 
as well as the asymptotic behaviour at large S (the tail of P{S)), upon the main 
features of the matrix element distribution functions (7^(0), (76^(61), (763(62). In the 
special case of an ensemble of random symmetric matrices A = A"*" we must have 
61 = 62. It is not enough that the two statistics gb-^ and gb^ are equal! This is 
achieved by inserting in the integrand of the constraint (7^2(62) = 5(6i — 62). 
Integrating then over 62 results in the level distribution formula 



PiS) 



da dbi 6 (s - 2\J a? + 6f J ga{a)gbiibi). 



(6) 



We shall work out exact formulae for the following cases: (i) General, non-normal 
matrices with Gaussian distributed elements 0,61,62 in sections |2] and [3], and (ii) 
symmetric (normal) matrices but considering non-Gaussian distributions of the ma- 
trix elements in section HI uniform (constant or box) distribution, Cauchy-Lorentz 
distribution, exponential distribution, and singular distribution (integrable power 
law at S" = multiplied by an exponential tail). Here we also detail more on the 
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Gaussian case in order to extend our results of reference f7] . In section [5] we com- 
ment on the level distribution of the prototype non-normal matrix, the triangular 
matrix. The final section [H] is devoted to a discussion and conclusions. 



2.2 Polar coordinate representation of the level spacing dis- 
tribution 

We notice that the radicands in the arguments of the delta functions in both equa- 
tions ([5]) and ([6]) are homogeneous in the moduli of a, bi, and 62- Thus it is natural 
to introduce spherical or plane polar coordinates. In the general case, described by 
equation (I5l), we define 

bi = r cos6 cosip, b2 = r cos6 simp, a = rsm6, (7) 

where r G [0, cxd), —71/ 2 < 9 < 7r/2, and < < 2tx. Then we get for the level 
distance 



S 



a2 + 6162 



2r Q with Q{d,^) 



I sin^ 9 -\ — cos2 9 sin 2lp 



The Jacobian of the coordinate transformation is cos 9 and therefore da dbi db2 - 
r'^dr cos 9 d9 d(p. The r-integration can be carried out in favour of S resulting in 



_ Z"^" cos9 d9 dip 



-tt/2 Jo 

s s s 

Now, if the value of the double integral were regular at 5 = 0, the level repulsion 
would be quadratic, P{S) oc S*^. But Q has zeros at ^ = together with ip = 
0, 71/2, vr, 37r/2, 2tt, which are scanned in the integration. That leads to singularities 
of the integrand, explicit ones (l/Q^) as well as implicit ones (in the gf-arguments). 
In particular the in the g^s enforces a scan of the full matrix element distribution 
functions including their tails, for any nonzero value of S. One cannot simply expand 
the g^s in terms of S and so find the small S behaviour. Thus the picture of the 
^'-dependence is far from simple. The small S behaviour depends on all details of 
the matrix element distribution functions g including their tails. Therefore in this 
general case of independent a, bi, 62 we shall choose another approach to attack this 
problem and analyse it at least in the case of Gaussian g^s in section |3l 

In the case of real symmetric 2D matrices the transformation to plane polar coordi- 
nates is simpler and thus much more useful. Here the r-integration does not produce 
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singularities, since the analog of Q here is just 1. Introduce plane polar coordinates 
into equation ([H]), 

a = rcos(y9, hi = r simp, (10) 

where r G [0, oo) and < < 2ti. The Jacobian is r, i. e., da dbi = rdr dip. The 

level distance reduces to 2 \J a? + 1)1 = 2 r, which is independent of ip> in contrast to 
the 6*, ip dependent factor Q in the general case of equation ([H]). We now can do the 
r-integration immediately and get 

S r'^'^ S S 
^(•5) = - dip gai- cos ip) gb,i- simp). (11) 

Here, for S ^ one can use the power law expansions of the matrix element 
distribution functions ga,bi in terms of S, independent of the nature of the (?'s, 
Gaussian or non-Gaussian. In particular, if ga{x) and gtiix) are both regular at 
X = 0, the integrand at = is just a constant and equal to 5'a(0)5'bi(0). We obtain 
for small S 

P{S)^S-^ga{0)gbAO). (12) 

Thus in case of regular g^s at zero value we always have linear level repulsion P{S) oc 
S for real, symmetric, random matrices, whose amplitude reflects the probability 
density ga,bi{0) to find the matrix elements a = and 6i = in the matrix A. Higher 
order corrections in S can be derived by Taylor expanding the g^s around x = 0. 
We conclude that regular matrix element distribution functions ga^^i transform into 
regular level spacing distributions P{S). From equations ffTTl) and ffT^ we also notice 
that the level repulsion is not linear if ga{0) and gtiiO) do not exist, i. e., if the 
distributions ga{x) and gbi{x) are singular at x = 0, if there is infinite probability 
density for the matrix elements to have values x = a = OoTX = bi = 0. We shall 
study this important case in section HI 



3 General non-normal real 2D matrix ensemble 
with Gaussian distributed matrix elements 

3.1 Level spacing distribution P{S): general 

We start with equation ([5]) and observe that the dependence of the level distance 
S in the integrand of P{S) on bi and 62 is only through the product B = 6162- 
Therefore it is natural to introduce hyperbolic coordinates defined as 
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B = 6162, = 775 equivalent to = Bv, h\ = B/v. (13) 

Both B and v run over the entire interval (— cxo,oo), but always have the same 
sign, sgn B = sgn bi ■ sgn 62 = sgn v. Positive B or v indicate that both non- 
diagonal elements are positive or that both are negative. Negative values of the 
variables B and v, instead, describe the case of non-diagonal elements with different 
sign. The Jacobian determinant is J = l/(2|t'|), and for the area elements we have 
dbi db2 = dB dv/{2\v\). 

To analyse the integral further we assume Gaussian distributed matrix elements, 
though with possibly different variances cr^, af, and cig, for a, 61, and 62, respectively. 



ga{a) = — ^exp( -), gh^k) = ^exp( — |), i = l,2. (14) 

All three distributions are normalized to one. With this assumption the equation 
for P{S) obtains the form 



PiS) 



1 



00 roo rco da dB dv 



Bv B 



2\v\ 



exp - + + 



Va'^ + B 



(15) 



The integrand is even in the variable a, we thus can use da ^ 2 da. Next, 
the level distance delta function does not depend on v explicitly. But since the 
signs of V and B are coupled, only the first and the third quadrant of the {B,v)- 
plane contribute. In both these (-B, f )-quadrants it is bf = Bv = \B\\v\ > and 
62 = B/v = \B\/\v\ > 0, guaranteeing the convergence of the (5, w)-integrals. In 
the first quadrant we have positive B, while in the third one — is relevant for the 
level distance. This leads to 



PiS) 



/■°° , f a'^\ f°° dv f Bv 

/ da exp / dB — exp ^ 

Jo \ ^ / "'0 ^0 2v \ af 



Bj 



bis- 2Va2 + b\ + (5 f 5 - 2 ^/~ 



B 



(16) 



Now all variables a,B,v have to be integrated over positive values only. The sum 
of the delta-functions then is independent of the variable v and the f-integral can 
be performed, cf. [H] No. 3.478,4. 



~ B B 1\ ^/2B\ 



2v 



a 



an V 



(17) 
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Here Ko{x) is the modified Bessel function of second kind and zero order. Its 
argument is always positive, B > 0. The level spacing distribution P{S) becomes 



oo 



P{S) = ^/ dB Ko{^)-G{B), 



/ 2B 

I — 3 / """" ^^\^\ 

where G{B^ is the Gaussian averaged level distribution for fixed, given product 
B > of the nondiagonal elements but varying diagonal elements a. 



G{B) 



bis- 2Va2 + 5) + (5 (S - 2 Va^ - 5 ) , with5 > 0. (19) 




a"' 



a2 



X 

The calculation of GiJ^) is easy because of the delta function and can be done 
analytically. Consider the first delta-function. Its argument, denoted as /(a) = 
S — 2^/a? + B, with positive B, corresponds to real eigenvalues of the matrix A. 
The zeros Oj of the delta function contribute to the integral only if they are real 
and positive. There is only one real, positive Oj, provided B < 5*^/4. Then the 
variable u = AB / S"^ fulfils < n < 1 and the /-zero reads Oj = ^ ^- The 
weight of the delta function contribution is given by the inverse of the derivative of 
/, which is \f\ai)\~^ = {2\/l — u )"^. The first delta-function in G{B) then leads to 
G+re(-B) = exp(— ^(1 — m))/(2a/1 — u ). The label {+re) indicates that the term 
describes the case +B and real (re) eigenvalues. Transforming the variable B ^ u 
this G+re{B) contributes the following term to the level level spacing distribution 



P.re{S) = f K, (^) exp f-^(l - .)) . (20) 

This integral has been considered already in ref. [7j and leads to a level repulsion 
exponent p = 2 — Oiog, see also section 13. 2[ It will turn out that this contribution 
( l20l) is subdominant relative to the other two integrals for P{S), which in contrast 
to ( l20l) will lead to the repulsion exponent p = 1, cf. again section [321 

We now calculate the G(i?)-contributions coming from the second delta function. 
They are labelled by a (— ) sign (since —B enters) and correspond to real as well 
as to imaginary eigenvalues of the matrix A, depending on the size of B. They are 
labelled therefore by (re) if — 5 > and by {im) if — B < 0. The relevant 
(positive) zeros Oj of the argument of the delta function are = |^ Vl + u for all 
positive B or u = AB / S"^ > in case of (re) and = ^^/u — 1 for all positive B 
with B > 5*^/4 and thus all u = AB/S"^ > 1 in the case {im). The weights \f'{ai)\~^ 
of the delta function contributions are obtained from |/'(aj)| = 2 a/1 -|- u for the case 
(re) and from |/'(aj)| = 2a/m — 1 for [im). These formulae lead to the following two 
contributions to the level spacing distribution P{S): 
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and 



AS) 



du 



Kn 



2(71 (72 



exp 



52 



4(t2 



[u-l] 



(22) 



These three integrals (pUj) . and (1^ can be summed up to give the complete 
level spacing distrribution P{S). To do this one introduces variables y such that 
in all three cases the exponential is exp(— y) with y > 0. In the case {—im) one 
in addition substitutes y —y. The real eigenvalues, represented by equations 
fl^Uj) and fl^Tj) . give y-integrations from to S'^/Aa'^ and from S'^/Aa'^ to 00. The 
imaginary eigenvalues lead to an integral from to cxo or, equivalently, from — cxd to 
0. Respecting the always positive argument of Kq one evaluates for the complete 
level spacing distribution function in closed form 



P{S) 



s 



dy 



01O2 



4(t2 



exp[-\y\ 



(23) 



This formula has been reported to us independently by Professor H.-J. Sommers 
jl8j . Because the integral converges also for = 0, we can conclude from this 
formula that in leading order in S we have linear level repulsion due to the explicit 
factor S, namely 



PiS) = S ■ ^ r ^ Ko (^y) eM-y) + h-o-t. (24) 

The S"- independent integral can be expressed in terms of the hypergeometric function 
F as will be shown in equation ( l33l) of section 13.21 

From equation (1231) one might wish to work out also the higher order terms in S, 
stemming from the integrand. However, if at small S one formally expands the factor 
Kq in fl2^ in terms of an S'^-series, one obtains contributions to the integrand, which 
are no longer integrable. Therefore the small S behaviour of the y- or -u-integrals is 
highly nontrivial. This agrees with our earlier observation in section 12.21 We have 
to analyse that in detail by studying the individual integrals P+re{S), P-re{S), and 
P-im{S) of equations (^U^, and (1^ . respectively. 
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3.2 Level spacing distribution P{S): Details of small S be- 
haviour 



As we have seen in sections 12.21 and 13.11 for non-symmetric (and thus non-normal) 
matrices, the behaviour of P{S) at small 5* is very delicate and depends on the 
details of the distribution functions gafiX^) for the matrix elements. In order to 
achieve understanding we go back to equation ( l20l) and make the substitutions 
1 — M — i> m' — > M and S'^u/ (2(Ji(J2) = u' ^ u yielding 

C re fir,. 

P^reiS) = ^ J_ . / ^ Me - u) e-^\ (25) 
Here we have introduced the notations 



Z(TiCT2 Zcr 

If e is very small, e ^ 1, we can use the approximation of Ko{z) at small z, which is 
(dlj. No. 8.447,1 and 3, and 8.362,3), Koiz) = -lnf-C + 0(z), where C is Euler's 
constant 0.577215 . . .. Then the leading term in the above integral, also taking into 
account the Taylor expansion of e~"^", can be written as 



/■^ nil 

hre{e)= / -=(-ln(e-^)), (27) 
Jo \Ju 

which after a simple substitution t — u = v!^u can be found in [17] (No. 2.727,5). 
After the evaluation for small e we get 

/+,e(e) ~ 2v^(- In e) oc 5 In S''^ , (28) 
meaning that including the explicit factor S in (1251) we have 



P+,e(^) oc ^^InS-^ (29) 

The level repulsion exponent p for the case of positive B (which also guarantees 
real eigenvalues of the matrix A) therefore is p = 2 — O^o^. Positive B means 
that we consider Gaussian distributed non-normal, real 2D matrices, which have 
only positive or only negative non-diagonal matrix elements 61,62- It is under this 
constraint that they have the rather strong repulsion exponent p = 2 — 0/og, as 
reported already in reference [7]. 

The other two contributions (1211) and (1221) are different in behaviour. Indeed each of 
them gives rise to weaker level repulsion, i. e., to a smaller exponent p. They both 
exhibit linear level repulsion P(S') oc 5, thus p = 1. This dominates the stronger, 
quadratic repulsion for small S valid for Pj^re{S), so that the total P{S) has linear 
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level repulsion P{S) oc S, as obtained in (l23l) and (IMll . The analytical reason for 
the quite different S'-dependence of P^re{S), from equation I^U^, in contrast to that 
of P_re and P-im, according to equations ([21]) and ([221), ^^e finite range of the 
B- or M-integration in ( !20l) versus the infinite integration intervals in the cases of 
negative B, namely ( 12T|) and ( l22l) . Because of these infinite integration intervals 
one cannot Taylor expand the ii'o-function in the cases of P-re and P-im- Instead, 
its complete functional form including its tails affect the convergence and thus the 
S'-behaviour of the integrals in ([21]), ([22])- 

In order to demonstrate this explicitly, we use the same substitution as before and 
arrive from equation ( |2T1) at the following expression 



5exp(-^) foo du ,,..Au 



P-re{S) = —-^^^ ^==K,{u)e-^\ (30) 



with the relevant integral 



^ roc rj/^i 

I_reiS) = / —=Koiu)e~^\ (31) 
Jo \/u + e 

The meaning of e and A is the same as in equation fl2U|) . The integrand decays 
exponentially for large u since also Kq[u) does so. Thus the upper limit of the 
integral converges safely, independent of e. For small u the integrand can be es- 
timated in a similar way as before. The result is that the integral converges to a 
finite value for e ^ and not to zero, as in equations ( 1271) and ( l28i) for l^re[S\ i- 
e., I-re{S — » 0) 7^ is finite as S ^ 0. Consequently, from the explicit factor S 
in equation (IHUll . we conclude that the level repulsion is linear, p = 1. Note that 
in contrast to the finite integration range in the integral (|271) . which for small S 
behaves oc 5* log 5"^, the infinite range integral (J3T1) has a finite, nonzero limit for 
S ^0. 

It remains to estimate the contribution from the imaginary eigenvalues given by (l22l) . 
Again, similar substitutions of the integration variable u lead us to the expression 



Q roo fill 

P-rm{S) = ^ ^3 / —K,{t + u)e-^\ (32) 



e and A are defined in equations (1261) . As before, the integral has a finite value, 
which does not go to zero with e oc ^ 0, and thus the level repulsion of the 
contribution P-im according to equations (l22l) and ( |32l) is again linear, p = 1. In 
fact, the integral in (1321) at e = can be calculated according to [T7] (6.621,3) as 



°°*A-„Me-^^ffii;i;4^). (33) 



/o ^ ^ ' yrTA V2'2' 'A + 1 

where F is the hypergeometric function. 
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The conclusion is that the level repulsion of the complete level spacing distribution 
function ^(5*) is linear, p = 1, due to the contributions P^re{S) and P_jm(5'), 
whereas the contribution P^re{S) alone has the level repulsion exponent p = 2 — 0;og, 
as we have already found in [7J for the case of Gaussian random matrices with only 
positive or only negative nondiagonal elements. Note that the criterion for the 
different small S behaviour is, if the product 6162 = -B is always positive, i? > 0, or 
if there are also -B < contributions. The eigenvalues according to this latter case 
are apparently more abundant and thus have less repulsion (oc S), while the former 
ones with only a positive product 6162 lead to the repulsion behaviour oc S'^logS'"^. 

4 Symmetric (normal) real 2D matrix ensembles 
with various distributions of the matrix ele- 
ments 

In this section we treat 2D real random matrices which are symmetric, i. e., 
hi = 1)2 = h and thus i? = 6^ > always (from here onwards we drop the la- 
bels 1 and 2 in fei = 62)- Such matrices are normal, i. e., the matrix commutates 
with its adjoint. The generalization here is, that we allow for a broad variety of 
matrix element distribution functions ga,b- The following classes of matrix element 
distributions are considered: (1) Gaussian distribution revisited, (2) box (uniform) 
distribution, (3) Cauchy-Lorentz distribution, (4) exponential distribution, and (5) 
singular distribution (power law approaching zero value, multiplied by exponen- 
tial tail). In all these cases we shall start with equation (ITTl) as a useful integral 
representation for ^(5*) in terms of Qa^- 

4.1 Gaussian distribution revisited 

This case of real, Gaussian distributed 2D matrices with possibly different widths 
of the diagonal and the non-diagonal matrix element statistics has been treated 
recently in reference [7]. We briefly repeat it here for the sake of completeness. 
Using the normalized Gaussians defined in (fT^ we immediately get 

where Iq{z) is the modified Bessel function of the first kind and zero order. According 
to ^J(No. 8.447,1) its small-z expansion is Iq{z) = 1 + x + fi + 0{z^). Thus 
/o(0) = 1 and the level repulsion is linear, p = 1. The details for the general case of 
different widths of the a and b statistics, cTq 7^ ab, has been analyzed and discussed in 
[7] . We mention that in the case of equal statistics of all matrix elements aa = (Jb = a 
we get, of course, the well known 2D GOE result 
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(35) 



After normalizing the first moment to unity, (S) = 1, leading to a = I/a/vt, 
the level spacing distribution P{S) becomes the Wigner distribution Pwigner{S) = 
^exp(-7r5V4). 

4.2 Box (uniform) distribution 

Let us now consider very different matrix element distributions ga,b- We start by 
studying the following uniform distributions for the matrix elements a and b, 

9a{ci) = - — , if \o\ < Go, otherwise , (36) 

gb{h) = —J-, if \b\ < bo, otherwise . (37) 

Thus the probability density product ga{,ci)gh{b) is constant, equal to (4ao&o)~^5 inside 
the (centrally located) rectangle with sides 2ao x 26o, and therefore for 5* smaller 
than 2min{ao,&o} the level spacing distribution P{S) can be calculated exactly. It 
is equal to 

nS) = ^. if ^<2min{ao,6o}. (38) 
8aoOo 

Again the level repulsion is linear. This is consistent with our finding in section 12.21 
that whenever ga,b{^) = const ^ there is generically linear level repulsion p = 1. 

From this geometrical picture it is also clear that P{S) is zero for S > 2^JaQ + feg- 
For S in between P{S) varies continuously. 

Indeed, we can calculate P{S) exactly for all S, using equation (fTT]) . We observe that 
ga and gb enter this expression symmetrically. Therefore without loss of generality 
we assume that bo < Oq, i.e. min{ao, 60} = ^o- Then we have to consider two more 
intervals, namely (i) 2bo < S < 2ao and (ii) 2ao < S < 2\JaQ + 6q. 

In the first case (i) the circle of radius S/2 intersects the rectangle at four points. The 
angle enclosed by the polar ray and the abscissa is equal to (po, where S sin ipo = 26o- 
Therefore, the total length of the 99-interval contributing to the integral in ffTTl) is 
just 4(y9o, and consequently P{S) = ^^|^arcsin^. 

In the second case (ii) the circle of radius S/2 intersects the rectangle at eight 
points, each pair of intersection points defining one y9-interval where we get the 
contribution to the integral. But all four angles are the same due to the double 
reflection symmetry of the rectangle and circle. The larger angle (p2 between the 
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polar ray and the abscissa is geometrically determined by Ssm(p2 = 26o, and the 
smaller one by Scos^pi = 2aQ. Thus for each pair the length of the yj-interval is 
equal to — ipi, and since there are four such intervals, the total length of the 
interval contributing to the integral is 4(y92 — V^i)- 

Putting all together we obtain the following exact result for the level spacing distri- 
bution P{S) in the case of the uniform (box) distributions ga,b, 



P{S) 



as < 2bo < 2ao, 
arcsin%^, if 2feo < 5 < 2ao, 



Sanbo 



S 



(arcsin ^ - arccos ^) , if 2ao < 5 < 2yagT&g, 



(39) 



4ao6o 

0, ifS> 2Jal + bl 



If instead of Bq < one has bo > oq, simply interchange oq and bo in the above 
formulae. 



4.3 Cauchy-Lorentz distribution 

The normalized probability densities for the matrix elements are defined by 

9a{a) = \ ^2. , ga{b) = ^ , . (40) 

"o "o 

From equation (fTTll and using plane polar coordinates we obtain 



P(s) = ^ r ^ (41) 

An^aoboJo (l + ^cosV)(l + ^sinV) ■ 

The integral for ^ gives 27r, so that at small S we have P{S) ~ S'/(27rao&o) in 
accordance with equation (fT2l) . But the integral (HTl) can also be done exactly. The 
result is 



P(S) = ^- «^yi+^ + /^^Vl+^ .42) 

^ ^ 27rao&o (a^ + + a^P^)VTT^VTTp^ ' 

where = S'^/^Aal) and = 6*^/(460). The asymptotic behaviour of P{S) at 
large S, i. e., ^ 1 and ^ 1, is an inverse quadratic power law, 

PiS) « l^f^. (43) 
If ao = bo = a, the complete formula for all level distances S reads 
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P(S) = , ^ , with = SVUa^). (44) 

This expression for the level spacing statistics evidently mirror images the Cauchy- 
Lorentz distribution of the {gt = (7a)-statistics in the P(S')-statistics. 

It is interesting to note that P{S) in has a divergent (infinite) first moment, as 
is clearly seen from the asymptotics fHHj) . A generalized power law statistics of the 
type gaia) = Ca/ (1 + (a/aoY), with g = 4, 6, . . ., however, has a finite first moment 
(5*) < oo. It will be treated elsewhere. 



4.4 Singular times exponential distribution 

Our normalized distributions in this subsection are chosen as 

ga{a) = a|a|-^''e-^'"l'^l, gb{b) = ai^r^'e-^^l^ (45) 
where the normalization constants are 



a = Al-^7(2r(l -/.,)). (46) 

Here i = a,b, the exponents /ij < 1, and T{x) is the gamma function. These 
distribution functions are singular but integrable power laws for a, 6 — and decay 
nearly exponentially in the tails. Using equation fill I) and the reflection symmetry 
(evenness) of both distributions ga^b in (H^ we get (note that S is positive only, 
^ > 0) 

P{S) = CaaS[-j cosM. ^ sin'^^ ^ • ^^^^ 

We did not succeed to calculate this integral analytically in closed form. However, 
one can evaluate it for small argument S, i. e., S ^ 0, where the exponential can 
be approximated by 1 (equivalent to A^, A^ 0, i. e., no tail effects in this small S 
range). The integral then is 



dip _ r(|-/i,)r(i-/ife) 



(4^ 



cos'^- if sin'"' if 2r(l - ^ - ^) 

From this we get the following level repulsion law, now being a fractional exponent 
power law, 

P(S) = 5 (1)-'"""" ^'4" "^'^'-I'^'l^' . (49) 
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The power law distribution for the matrix elements leads to a corresponding power 
law for the level spacing distribution. The power law exponents and fit of the 
matrix element distribution functions Qa^b immediately transform into the level re- 
pulsion exponent. More precisely, the level repulsion exponent is p = 1 — /Iq — /Xb. 
We emphasize that P{S) at S" = is integrable, if the matrix element distributions 
ga{ci) and gb{b) are integrable at a = and 6 = 0, respectively. 

The physical interpretation of this repulsion exponent p = 1— fj,a — Pb comprises two 
different possibilities, depending on the size of the singularity exponents and Hb 
of the distribution functions Qa, Qb for the matrix elements. If the singularities are 
strong, more precisely, if /ia + /ib > 1, the repulsion exponent p is negative. This 
means that due to the rather strong sparsing of the matrix there is not a repulsion 
but, instead, an enhancement of the level distance S = 0, the zero eigenvalues 
are emphasized. If the singularities are weak, pa + f^b < 1, the phenomenon of 
level repulsion remains, p > 0, despite the singularities for the matrix element 
distributions at a = and 6 = 0. 

The diagonal elements ±a and the non-diagonal elements b determine the level 
distance with equal weight since S = 2y/a?+l^. It is for this reason that it is just 
the sum of the singularity exponents which determines p, giving equal weight to 
the singularities of the diagonal and non-diagonal elements to the repulsion. There 
are two typical limiting cases, (i) pa = Pb = p, both singularities are of equal 
strenght, or (ii) pa = po ^ and pb = or vice versa. Only one of the two matrix 
element distributions is singular while the other one remains regular. Then only 
the non-diagonal is sparsed and the diagonal is regular or the other way round. In 
the first case (i), if /i < 1/2 we have still level repulsion, p > 0, while for p > 1/2 
(but still less than 1) we find enhancement at zero distance 5 = 0. - In the second 
case (ii) there will always be level repulsion despite the singularity at zero for either 
the non-diagonal or the diagonal distribution, since p = 1 — po > always. The 
repulsion exponent in this case will be between and 1. Sparsing like this is one of 
the possible causes that the repulsion exponent scans the interval (0, 1) . 

These results are interesting in the context of quantum chaos of nearly integrable 
(KAM type) systems. As has been observed in a variety of different systems of 
mixed type, at small energies one finds the so-called fractional exponent power law 
level repulsion, well described by the Brody distribution rather than by the Berry- 
Robnik distribution ^19j-[2l], which in turn has been clearly demonstrated to apply 
at sufficiently large energies, i. e., at sufficiently small effective Planck constant. The 
observed deviation from the Berry-Robnik behaviour is due to the localization and 
tunneling effects, and is a subject of intense current research [25]. Phenomenologi- 
cally it has been discovered in [8], further developed in [9], [lO], and the connection 
with sparsed matrices was established in pjj. Quite generally, the matrix repre- 
sentation of a Hamilton operator of a nearly integrable (KAM type) system in the 
basis of the integrable part results in a sparsed banded matrix with nonzero di- 
agonal elements [TT]. Such a sparsed matrix is precisely characterised by the fact 
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that many matrix elements are zero. In other words, the probabihty distribution 
function of the nondiagonal matrix elements gb{b) is singular at zero value of b, in 
a manner described by equation fH31) . while the diagonal matrix elements have a 
regular distribution function, which is precisely the case (ii) discussed above. Our 
2D random matrix theory with such singular matrix element distribution functions 
therefore predicts qualitatively a fractional exponent power law level repulsion, the 
phenomenon observed in the above mentioned works [8]- [11]. Thus we see that the 
study of random matrices with other than the invariant ensembles (GOE and GUE) 
is very important and connects to new physics. We leave this direction of research 
for further studies in the near future. 

The large S behaviour of P{S) in this case is obviously dominated by the expo- 
nentials. Although an exact solution of the integral ( 1471) is not known, we shall 
show by applying the mean value theorem to the integral that P{S) decays roughly 
exponentially at large S ^ 1. This will be analysed in the next subsection, devoted 
to pure exponential matrix element distributions without power law singularities at 
the origin. 



4.5 Exponential distribution 

Here we start with the distribution functions of the previous subsection, but without 
singularities. I. e., we assume fJ-a = fJ'b = 0, yielding the case of purely exponential 
distribution of matrix elements. 

ga{a) = ae-^'^l'^l, gb{b) = ae"^"''', with = X,/2. (50) 
The level spacing distribution function in this case is 

P{S)=CaCbSj^ dip exp i- — {Xa COS (f + Xb sin (f)] . (51) 

We could not evaluate this analytically in closed form. For small S the linear level 
repulsion law with p = 1 is recovered, of course, 

PiS) ^ ^ S=^S. (52) 
2 8 

At large 5* we can estimate the integral fl^T]) using the mean value theorem. First we 



write Xa cos (f + Xb simp = A sin{ip + (p) , with = arctan(Aa/A;,), and A = yA^ + Af. 
We then substitute the integration variable from <^ to x = </9 + 0, which now runs 
from to + 7r/2. This transforms P{S) into 

P{S) = CaCb S / dx e-^^'"'^. (53) 
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Since the integrand is continuous and bounded we can apply the mean value theorem, 
saying that there is a value Xo('S') in the interval between and + 7r/2 such that 
the level spacing distribution is 

P(5) = . s ■ e-P--xo{s)_ (54) 

Xo{S) is expected to be only weakly dependent on S. In this sense the tail of P{S) 
is roughly exponential. 



5 Triangular matrix 

A special, prototype case of a non- normal matrix A from equation ([T]) is the trian- 
gular matrix 




(55) 



Such matrices have been studied in detail in reference [7]. For them 6162 = B = 
and the level distance S does not depend on b. With the constraint (762(^2) = <5(^2) in 
equation ([5]) the a and (61 = 6)-integrations factorize, leading to the level distribution 
function 



P{S) = / da 6{S -2\a\)gaia) 
For even distributions ga{ci) we have 



S 

ga\- 



s 

9a\-- 



(56) 



P{S) = ga 



S >0, normalized on [0 < 5 < 00) 



(57) 



Apparently for triangular matrices the level spacing distribution is the immediate 
mirror of the diagonal element distribution ga{a)- If ga{o-) is Gaussian, there is no 
level repulsion at all, the repulsion exponent is p = (see already [7J), though still 
the S are Gaussian and not Poissonian distributed as in the integrable case. If ga{ci) 
is Poissonian, so is P{S). And if ga{ci) is centered about some finite value a, this 
also holds for the level spacing distribution. If the a-distribution does not include 
the origin, there even is formally infinite level repulsion. For general matrices with 
non-zero mean values d and 61^2 the level spacing distribution is somewhat more 
tricky. 
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6 Discussion and conclusions 



We have studied the level spacing distribution of general 2D real random matri- 
ces. First general, non-normal random matrices with Gaussian distribution of the 
matrix elements have been considered, showing that in general the level repulsion 
is still linear, as for symmetric matrices. Under some constraints it is quadratic 
with logarithmic corrections. We have given an explicit formula for the level spac- 
ing distribution function in form of a threefold integral. Then we have considered 
symmetric matrices with general, other than Gaussian statistics for both the diag- 
onal and the non-diagonal elements. We have shown that the level repulsion again 
is always linear provided the matrix element distribution functions are regular at 
zero value with finite, non-zero weight. We have explicitly considered the box-type 
(uniform), the Cauchy-Lorentz, the exponential times singular power law at zero, 
and the purely exponential matrix element distributions. Explicit closed form re- 
sults for P{S) have been obtained, except for the singular times exponential tail 
distributions, although we could present a very good understanding of the overall 
behaviour of P{S) also in this case. 

Our approach can obviously be extended to general 2D complex random matrices, 
first to Hermitian complex matrices, where similar general formulae like equation 
( fTTj) can be obtained, by using spherical coordinates. Indeed, it becomes obvious 
that we shall always have quadratic level repulsion as long as the matrix element 
distributions g{x) are regular at zero x, and a formula analogous to ( fT2l) will ap- 
ply. Denoting the diagonal elements by ±a, and the off-diagonal elements (complex 
conjugate) by 6 ± z c, we obtain 



S"^ f f S S S 

P{S) = — / cos 6 d9 dip ga{— sin 6) gb{— cos 9 sin (p) gc{— cos 6 cos (f), (58) 
8 J J 2 2 2 

as claimed. The details of further analysis will be published in a separate paper. 

We have also discussed the connection between the singular distributions of the 
matrix elements with a power law at zero value and the sparsed matrices which 
describe nearly integrable systems. In such systems one derives a fractional exponent 
power law level repulsion well known phenomenologically but poorly understood 
theoretically. Our approach and results promise new advances in this important 
direction of research in quantum chaos of mixed type systems initiated in [2S] and 
general random matrix theory. 
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